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Abstract 

Functional interpretations of limb bone structure frequently assume that diaphyses adjust their shape by adding bone 
primarily across the plane in which they are habitually loaded in order to minimize loading-induced strains. Here, to test this 
hypothesis, we characterize the in vivo strain environment of the sheep tibial midshaft during treadmill exercise and 
examine whether this activity promotes bone formation disproportionately in the direction of loading in diaphyseal regions 
that experience the highest strains. It is shown that during treadmill exercise, sheep tibiae were bent in an anteroposterior 
direction, generating maximal tensile and compressive strains on the anterior and posterior shaft surfaces, respectively. 
Exercise led to significantly increased periosteal bone formation; however, rather than being biased toward areas of 
maximal strains across the anteroposterior axis, exercise-related osteogenesis occurred primarily around the medial half of 
the shaft circumference, in both high and low strain regions. Overall, the results of this study demonstrate that loading- 
induced bone growth is not closely linked to local strain magnitude in every instance. Therefore, caution is necessary when 
bone shaft shape is used to infer functional loading history in the absence of in vivo data on how bones are loaded and how 
they actually respond to loading. 
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Introduction 

Functional analyses of limb bone structure, without prior 
knowledge of loading history, commonly assume that the 
distribution of bony material within diaphyseal cross sections is 
affected by the magnitude of local mechanical strains (deforma- 
tions) generated by habitual physical activities such as locomotion 
[1]. In this strategy, bone tissue's formative and resorptive activity 
is considered an optimization process that acts primarily to 
maximize resistance to loading-induced strains with a minimum of 
material. Functional loading is expected to have a predictable 
effect on diaphyseal shape by promoting bone formation 
disproportionately in areas of the shaft surface that are routinely 
subject to the highest levels of strain. Thus, bone shaft cross- 
sectional shapes that display structural reinforcement in a 
particular plane (e.g., elliptically shaped cross sections) are 
hypothesized to have experienced more frequent bending loads 
and higher strains across the axis about which the bulk of bone 
tissue is distributed [2-4]. 

However, studies documenting in vivo diaphyseal strains in a 
variety of animals have shown that limb bone shaft cross sections 
are not always reinforced in the planes in which they are habitually 
loaded [5-8] . For example, during quadrupedal walking, the ulna 
of a macaque monkey primarily experiences coronal bending, with 



the highest strains produced on the medial and lateral shaft 
surfaces; however, macaque ulnar diaphyses are shaped to better 
resist sagittal bending [7]. Some functional morphologists have 
downplayed the relevance of these results by pointing to other 
strain studies that found the correspondence between maximum 
strain distribution and maximum bending rigidity in diaphyses to 
improve during more vigorous activities [1]. For example, when 
macaques transition from walking to galloping, the greatest strains 
engendered in tibial shafts begin to approximate the direction of 
maximal diaphyseal bending rigidity [9], This suggests to some 
that diaphyseal shape indeed reflects patterns of functional 
loading, particularly activities involving vigorous loading [1]. 

In order to rigorously evaluate the relationship between limb 
bone loading history and diaphyseal shape, experiments are 
required that directly relate in vivo functional strains to patterns of 
bone growth. To date, however, few studies of this kind have been 
conducted. In the present study, we characterize the in vivo strain 
environment of the tibial midshaft in sheep during treadmill 
exercise and test whether bone formation induced by this activity 
occurs primarily in diaphyseal regions that experience the highest 
strains. Such a growth response would structurally reinforce the 
shaft across its bending axis and, likewise, would influence 
diaphyseal shape in a manner consistent with the model within 
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Figure 1. Tibial midshaft cross sections. The sedentary control sheep (subject 4) and the exercised sheep (subject 3) had similar body masses 
throughout the experiment. The calcein labels injected after the first week of the exercise treatment are visible as the light grey rings within the 
sections. Scale equals 1 mm. 
doi:10.1371/journal.pone.0099108.g001 



which functional morphologists often interpret shaft structure in 
the absence of information on loading history. 

Materials and Methods 

Ethics Statement 

All experimental procedures were reviewed and approved by 
the Institutional Animal Care and Use Committee of Harvard 
University. 

Experimental Procedures 

Juvenile animals were employed because the osteogenic 
response to exercise is typically greatest during the growth period 
[10-12]. At 40 days of age, rams (Ovis dries; Dorset) were divided 
into three groups. In the first group (N = 5), strain gauges were 
attached to the tibial midshaft, and strain data were recorded 
during treadmill trotting at a Froude speed of 0.5 (approximately 
4 km h ). A second group (N = 5) was exercised on a treadmill at 
the same Froude speed for 30 min day -1 (approximately 6,000 
loading cycles day ) for 90 days. Subjects were injected (i.p.) with 
calcein (20 mg kg -1 ) after the first week of the exercise treatment 
to record bone growth on the periosteal surface of the tibia. The 
third group (N = 5) served as sedentary controls and received a 
calcein injection but no exercise treatment. All subjects were 
individually housed and were given the same ad libitum diet of food 
and water. After the 90-day treatment period, subjects were 
euthanized and tibiae were extracted. 

Strain Gauge Recordings 

Following induction and maintenance of isoflurane general 
anaesthesia, triple-rosette strain gauges (Sokki Kenkyujo, Tokyo, 
Japan) were affixed to the anterior, medial, and posterior 
diaphyseal surfaces of the left tibia through small skin incisions 
using aseptic technique. At each gauge site, 5 mm 2 of periosteum 
was removed, the bone surface was degreased with 100% 
chloroform, and a gauge was glued to the prepared surface with 
methyl-2-cyano-acrylate. Care was taken to align one gauge 
element with the long axis of the bone. Gauge leads were passed 
extracutaneously beneath flexible bandages to the hip, where they 
were sutured to a belt worn around the subject's abdomen. Strain 
relief was provided by attaching the leads to dressing wrapped 
around the leg near the incision sites. 

Strain data were recorded 4 and 24 h after surgery, when 
subjects were moving with normal gaits and showed no signs of 



lameness or discomfort. During each recording session, gauges 
were connected with insulated wire to Vishay 2 120 A amplifiers 
(MicroMeasurements Inc., Raleigh, NC) to form one arm of a 
Wheatstone quarter-bridge circuit. Bridge excitation was 1 V. 
Voltage outputs were recorded on a Teac RD-145T DAT tape 
recorder (Teac Corp., Tokyo, Japan). Gauges were calibrated 
when subjects were stationary with their instrumented limb off the 
ground, and were periodically balanced during experiments to 
adjust for zero offsets. 

Strain Gauge Analyses 

Strain data sequences were sampled from tape recordings using 
an Ionet A-D board (GW Instruments, Somerville, MA) at 
250 Hz. A custom-written Superscope 3.0 (GW Instruments, 
Somerville, MA) virtual instrument was used to determine the zero 
offset, calculate strains of principal tension and compression from 
raw voltage data using shunt calibration signals recorded during 
experiments, and calculate the orientation of principal tensile 
strain relative to the bone's long axis using standard engineering 
equations [13]. Igor Pro 4.0 (Wavemetrics Inc., Lake Oswego, 
OR) was used to calculate these strains at temporal midstance — 
when peak strains occur [8] — for at least 10 gait cycles per subject. 
To characterize the average tibial midshaft strain environment 
among all subjects, a custom-written macro [8] for ImageJ (NIH, 
Bethesda, MD) was used to compute the mean location and 
orientation of the neutral axis, as well as mean gradients of cross- 
sectional normal strain, using linear beam theory [13]. 

Histological Analyses 

In the exercised and control subjects, bone formation on the 
tibial periosteal surface during the treatment period was quantified 
on midshaft transverse cross sections (100-|tm thickness) using an 
Olympus SZH-10 microscope (Olympus America, Melville, NY) 
with epifluorescence (Fig. 1). Using ImageJ, digitized cross sections 
were divided into 16 equal-angle sectors positioned about the 
experimentally determined neutral axis, with the axis orthogonal 
to the neutral axis projected through the area centroid [8]. Bone 
area added in each sector during the treatment period was 
measured from the calcein line to the periosteal surface. Bone 
areas were standardized by body mass ' , calculated as mean 
mass during the final 3 weeks of the experiment. Two-sample 
Wilcoxon tests were performed to evaluate differences in bone 
formation between exercised and control animals. 
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Figure 2. Bone strain and bone formation. (A) Mean distribution of peak longitudinal normal strain (u,s) across the tibial mid-diaphysis during 
treadmill exercise. The neutral axis (NA) is dashed with compressive and tensile strain isopleths plotted parallel to it. Maximal compressive strains 
(negative sign) were produced on the posterior shaft surface, and maximal tensile strains were produced on the anterior shaft surface. (B) Tibial 
midshaft cross section subdivided into 1 6 equal-angle sectors positioned about the neutral axis, with the axis orthogonal to the neutral axis projected 
through the area centroid. (C) Distribution of periosteal bone added in exercised animals and sedentary controls during the experimental period 
(means + SD). Asterisks indicate statistically significant (P<0.05) differences between exercised animals and controls as determined by two-sample 
Wilcoxon tests. Numbers indicate percent difference between the mean of the exercised animals relative to the control mean. 
doi:1 0.1 371 /journal. pone.0099108.g002 



Additional details of the experimental and analytical methods 
are provided elsewhere [8,11]. Raw bone strain and histological 
data are available in [8] and Dataset SI, respectively. 



Results and Discussion 

Treadmill exercise caused a highly non-uniform strain environ- 
ment in the sheep tibial midshaft (Fig. 2A). During peak loading, 
maximal compressive strains were generated on the posterior shaft 
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surface, and maximal tensile strains were produced on the anterior 
surface. The neutral axis was oriented within 10° of a mediolateral 
axis and shifted toward the anterior diaphyseal surface, indicating 
that bones were bent in an anteroposterior direction. Maximal 
compressive strains were approximately 40% higher than tensile 
strains, as expected in a loading regime that includes both bending 
and axial compression. Based on these data, if functional loading 
stimulates bone formation disproportionately in diaphyseal regions 
that experience the highest strains, as is often assumed in 
functional interpretations of bone shaft structure, then one would 
expect treadmill exercise to promote osteogenesis primarily on the 
anterior and posterior tibial shaft surfaces. 

However, in the sheep in this study, locations of exercise- 
induced bone formation were poorly linked to local strain 
magnitude during peak loading (Fig. 2B,C). Instead of being 
biased toward areas of maximal strains across the anteroposterior 
shaft axis, significant exercise-induced bone formation occurred 
primarily around the medial half of the diaphyseal circumference, 
in both high strain regions and direcdy adjacent to the neutral axis 
where strains were lowest. In nearly all areas of the lateral shaft 
surface, exercise failed to stimulate significant bone formation, 
including in multiple sectors that experienced maximal strains. 
Relative differences in bone added between exercised sheep and 
sedentary controls were greatest in an area of high tensile strains 
(sector 2) and smallest in the two lateral sectors adjacent to the 
neutral axis (sectors 12 and 13), yet relative group differences were 
very similar between other regions of maximal and minimal strains 
(e.g., sector 4 vs. 8). In sum, these results demonstrate that, 
although it might be assumed that the primary purpose of bone 
response to loading is to adjust its morphology to best resist 
maximal strains, this is not necessarily what is achieved in every 
instance. 

These findings accord well with those of previous experiments 
with bird models involving both exogenous limb loading [14] and 
weight-bearing exercise [6] that failed to find any significant 
correlation between bone formation stimulated by loading and 
local strain magnitude. In contrast, however, peak magnitude 
strains have been associated with sites of increased bone formation 
in experiments involving external loading of rodent limbs [15,16] 
and natural physiological limb loading in goats [17]. A notable 
difference between these previous studies is that in the bird 
experiments, subjects were skeletally mature and the surfaces of 
their bones were primarily quiescent, whereas subjects in the 
mammal experiments were growing and had bones with active 
surfaces. Bone's response to loading is known to be affected by age 
[10-12], and has been hypothesized to differ between distandy 
related taxa [1,18]. Therefore, the agreement between the results 
of bird studies and the present experiment involving growing 
mammals is important because it implies that the bone functional 
response observed in sheep was not due to their specific age, nor is 
it unusual for vertebrates in general. Also important is the disparity 
between our results and those of other mammal studies because it 
suggests that a universal relationship between bone growth and 
local strain magnitude probably does not exist. 

Previous in vivo strain studies demonstrating that, in many 
animals, limb bone shafts are structurally reinforced about the 
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